Search for: All records

We investigate the local proton energization at magnetic discontinuities/intermittent structures and the corresponding kinetic signatures in velocity phase space in Alfvénic (high cross helicity) and non-Alfvénic (low cross helicity) wind streams observed by Parker Solar Probe. By means of the partial variance of increments method, we find that the hottest proton populations are localized around compressible, coherent magnetic structures in both types of wind. Analysis of parallel and perpendicular temperature distributions suggest that the Alfvénic wind undergoes preferential enhancements ofT_∥at such structures, whereas the non-Alfvénic wind experiences preferentialT_⊥enhancements. Although proton beams are present in both types of wind, the proton velocity distribution function displays distinct features. Hot beams, i.e., beams with beam-to-core perpendicular temperatureT_⊥,b/T_⊥,cup to three times larger than the total distribution anisotropy, are found in the non-Alfvénic wind, whereas colder beams are in the Alfvénic wind. Our data analysis is complemented by 2.5D hybrid simulations in different geometrical setups, which support the idea that proton beams in Alfvénic and non-Alfvénic wind have different kinetic properties and different origins. The development of a perpendicular nonlinear cascade, favored in balanced turbulence, allows a preferential relative enhancement of the perpendicular plasma temperature and the formation of hot beams. Cold field-aligned beams are instead favored by Alfvén wave steepening. Non-Maxwellian distribution functions are found near discontinuities and intermittent structures, pointing to the fact that the nonlinear formation of small-scale structures is intrinsically related to the development of highly nonthermal features in collisionless plasmas. Our results contribute to understanding the role of different coherent structures in proton energization and their implication in collisionless energy dissipation processes in space plasmas.

 

Revealing the formation, dynamics, and contribution to plasma heating of magnetic field fluctuations in the solar wind is an important task for heliospheric physics and for a general plasma turbulence theory. Spacecraft observations in the solar wind are limited to spatially localized measurements, so that the evolution of fluctuation properties with solar wind propagation is mostly studied via statistical analyses of data sets collected by different spacecraft at various radial distances from the Sun. In this study we investigate the evolution of turbulence in the Earth’s magnetosheath, a plasma system sharing many properties with the solar wind. The near-Earth space environment is being explored by multiple spacecraft missions, which may allow us to trace the evolution of magnetosheath fluctuations with simultaneous measurements at different distances from their origin, the Earth’s bow shock. We compare ARTEMIS and Magnetospheric Multiscale (MMS) Mission measurements in the Earth magnetosheath and Parker Solar Probe measurements of the solar wind at different radial distances. The comparison is supported by three numerical simulations of the magnetosheath magnetic and plasma fluctuations: global hybrid simulation resolving ion kinetic and including effects of Earth’s dipole field and realistic bow shock, hybrid and Hall-MHD simulations in expanding boxes that mimic the magnetosheath volume expansion with the radial distance from the dayside bow shock. The comparison shows that the magnetosheath can be considered as a miniaturized version of the solar wind system with much stronger plasma thermal anisotropy and an almost equal amount of forward and backward propagating Alfvén waves. Thus, many processes, such as turbulence development and kinetic instability contributions to plasma heating, occurring on slow timescales and over large distances in the solar wind, occur more rapidly in the magnetosheath and can be investigated in detail by multiple near-Earth spacecraft.

 

Abstract Several generalizations of the well-known fluid model of Braginskii (1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multispecies plasmas with arbitrary masses and temperatures, where all of the fluid moments are described by their evolution equations. The 21-moment model contains two “heat flux vectors” (third- and fifth-order moments) and two “viscosity tensors” (second- and fourth-order moments). The Braginskii model is then obtained as a particular case of a one ion–electron plasma with similar temperatures, with decoupled heat fluxes and viscosity tensors expressed in a quasistatic approximation. We provide all of the numerical values of the Braginskii model in a fully analytic form (together with the fourth- and fifth-order moments). For multispecies plasmas, the model makes the calculation of the transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasistatic approximation that makes some Braginskii coefficients divergent in a weakly collisional regime. Importantly, we show that the heat fluxes and viscosity tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the fourth-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad. 

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(\unicode[STIX]{x1D701})$ , and the associated plasma response function $R(\unicode[STIX]{x1D701})=1+\unicode[STIX]{x1D701}Z(\unicode[STIX]{x1D701})=-Z^{\prime }(\unicode[STIX]{x1D701})/2$ . We then consider a one-dimensional (1-D) (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the fourth-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1-D closures at higher-order moments than the fourth order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3-D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the fourth-order moment level. In appendix A, we provide comprehensive tables with Padé approximants of $R(\unicode[STIX]{x1D701})$ up to the eighth-pole order, with many given in an analytic form. 

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modelling the turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non-Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the fluid model of Chew–Goldberger–Low (CGL) in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the fourth-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the ‘normal’ closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the ‘normal’ closure, and beyond the fourth-order moment, Landau fluid closures are required. 

We use 2.5D magnetohydrodynamic simulations to investigate the spectral signatures of the non-linear disruption of a tearing unstable current sheet via the generation of multiple secondary current sheets and magnetic islands. During the non-linear phase of tearing mode evolution, there develops a regime in which the magnetic energy density shows a spectrum with a power law close to B(k)2 ∼ k−0.8. Such an energy spectrum is found in correspondence of the neutral line, within the diffusion region of the primary current sheet, where energy is conveyed towards smaller scales via a ‘recursive’ process of fast tearing-type instabilities. Far from the neutral line, we find that magnetic energy spectra evolve towards slopes compatible with the ‘standard’ Kolmogorov spectrum. Starting from a self-similar description of the non-linear stage at the neutral line, we provide a model that predicts a reconnecting magnetic field energy spectrum scaling as k−4/5, in good agreement with numerical results. An extension of the predicted power law to generic current sheet profiles is also given and possible implications for turbulence phenomenology are discussed. These results provide a step forward to understand the ‘recursive’ generation of magnetic islands (plasmoids), which has been proposed as a possible explanation for the energy release during flares, but which, more in general, can have an impact on the subsequent turbulent evolution of unstable sheets that naturally form in the high Lundquist number and collisionless plasmas found in most of the astrophysical environments.